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Abstract. In this paper, we continue to construct stationary classical solutions of the 
incompressible Euler equation approximating singular stationary solutions of this equa- 
tion. This procedure now is carried out by constructing solutions to the following elliptic 
problem 

I ~e 2 Au = YJ?=i Xn+ ( u ~ 1 ~ ln §)+ - E"=i Xaj (s " ^ ln \ - u Y+> x e n > 
\u = 0, xedCl, 

where p > 1, f2 C M 2 is a simply connected bounded domain, Slf and CIJ are mutually 
disjoint subdomains of Q, q is a harmonic function. 

We showed that if SI is a simply-connected smooth domain, then for any given (In- 
stable critical point of Kirchhoff-Routh function W{x 1 , ■ ■ ■ , x^, x 1 , ■ ■ ■ , x n ) with k+ > 
(i = 1, • • • , m) and nj > (j = 1, • • • , rt), then there is a stationary classical solution 
approximating stationary m + n points vortex solution of incompressible Euler equations 
with total vorticity Y^iLi K i ~ Ej=i K J ■ 
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1. Introduction and main results 
The incompressible Euler equations 

\ t + (v • V)v = -VP, 

V-v = 0, 1 ' ' 

describe the evolution of the velocity v and the pressure P in an incompressible flow. In 
1R 2 , the vorticity of the flow is defined by u = V x v := d±V2 — c^fi, which satisfies the 
equation 

u t + v • Vw = 0. 

Suppose that u is known, then the velocity v can be recovered by Biot-Savart law as 
following: 

1 -x L 

v — u * 
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where x 1 - = (x 2 , —x\) if x = (xi, x 2 ). One special singular solutions of Euler equations is 
given by u = YlT=i K ^Xi(t), which is related 

EKi {x — Xj(t)) -1 " 
, 2^ \x-xAt)\ 2 ' 
i=i 1 v ' 1 

and the positions of the vortices X{ : R — > M 2 satisfy the following Kirchhoff law: 
where W is the so called Kirchhoff-Routh function defined by 

m 

2 27T SCj — X 



:i l 



In simply-connected bounded domain ft C M 2 , similar singular solutions also exist. Sup- 
pose that the normal component of v vanishes on <9ft, then the Kirchhoff-Routh function 
is 

^ m ^ m 

■ • • , SC m ) = - 2j KiKjGfa, Xj) + - 22 f^H(Xi, Xi), (1.2) 

where G is the Green function of —A on ft with Dirichlet boundary condition and H is 
its regular part (the Robin function). 

Let v n be the outward component of the velocity v on the boundary <9ft, then we see 
that J dn v n = due to the fact that V ■ v = 0. Suppose that v is the unique harmonic 
field whose normal component on the boundary <9ft is v n , then vo satisfies 

V • v = 0, in ft, 

V x v = 0, in ft, (1.3) 
n ■ v = v n , on 9ft. 

If ft is simply-connected, then v can be written v = (Vil ! o)~ L , where the stream function 
ipo is determined up to a constant by 

-A^o = 0, in ft, 

v„, on oil 



Or 

where denotes the tangential derivative on <9ft. The Kirchhoff-Routh function associ- 
ated to the vortex dynamics becomes(see Lin [24]) 

W(xi,--- ,X m ) = ^ ^ KjKjGjXj, Xj) + n ^ K?H(Xj, Xj) + ^ Kj^Xj) ■ (1.5) 
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For m clockwise vortices motion (corresponding to nf > 0) and n anti-clockwise vortices 
motion (corresponding to — k~ < 0), the Kirchhoff-Routh function associated to the vortex 
dynamics becomes 

^ m 1 n 

i,k=l,i^k j,l=l,j^l 
1 m i n 

+ ^ K t) 2 H(xl,xi) + -J2^j) 2 H(xJ,xj) 

i=l j=l 
m n m n 

-YY K t K j G ( x ti x j) + Y K tM x t) - Y K jM x j)- 

i=l j=l i=l j=l 

(1.6) 

It is known that critical points of the Kirchhoff-Routh function W give rise to stationary 
vortex points solutions of the Euler equations. As for the existence of critical points of W 
given by ( ll.2j) . we refer to [5]. 

Roughly speaking, there are two methods to construct stationary solutions of the Euler 
equation, which are the vorticity method and the stream-function method. The vorticity 
method was first established by Arnold and Khesin [3] and further developed by Burton 
[7] and Turkington [32] . 

The stream-function method consists in observing that if ip satisfies —Atp = f(i/j) for 
some function / G C 1 (M), then v = (Vip) ± and P = F(ip) — ^{Vipl 2 is a stationary 
solution to the Euler equations, where (V^) -1 := {■ijjr, — F(t) = f f(s)ds. Moreover, 
the velocity v is irrotational on the set where f(ip) = 0. 

Set q = —ipQ and u = ip — tp , then u satisfies the following boundary value problem 

-Au = f(u- q), xett, ^ 
u — 0, x € dVt. 

In addition, if we suppose that inf^ q > and f(t) = 0, t < 0, the vorticity set {x : f{jp) > 
0} is bounded away from the boundary. 

The motivation to study fll.7j) is to justify the weak formulation for point vortex solu- 
tions of the incompressible Euler equations by approximating these solutions with classical 
solutions. 

Marchioro and Pulvirenti [26] have approximated these solutions on finite time intervals 
by considering regularized initial data for the vorticity. On the other hand, the stationary 
point vortex solutions can also be approximated by stationary classical solutions. See e.g. 
[H El H El UM ESI ESI EH E21 E31 El] and the references therein. 

In [T5] Elcrat and Miller, by a rearrangements of functions, have studied steady, inviscid 
flows in two dimensions which have concentrated regions of vorticity. In particular, they 
studied such flows which " desingularize" a configuration of point vortices in stable equilib- 
rium with an irrotational flow, which generalized their earlier work for one vortex [16j[17j 
which in turn were based on results of Turkington [32]. As pointed by Elcrat and Miller, an 
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essential hypothesis in their existence proof was that the vorticity was in a neighborhood of 
a stable point vortex configuration. Saffman and Sheffield [30J have found an example of a 
steady flow in aerodynamics with a single point vortex which is stable for a certain range of 
the parameters. This has been generalized in [16], where some examples computationally 
of stable configurations of two point vortices were briefly discussed. Further examples of 
multiple point vortex configurations are given in [27] , where a theorem on the existence of 
such configurations is also given. 

It is worth pointing out that except [18] the above approximations can just give explana- 
tion for the formulation to single point vortex solutions. D. Smets and J. Van Schaftingen 
investigated the following problem 

= (u - q - ^-ln-) p , in Q, , s 

on dQ, 

and gave the exact asymptotic behavior and expansion of the least energy solution by 
estimating the upper bounds on the energy. The solutions for (11.81) in [31] were obtained 
by finding a minimizer of the corresponding functional in a suitable function space, which 
can only give approximation to a single point non- vanishing vortex. In [13], we have shown 
that multi-point vortex solutions can be approximated by stationary classical solutions. 

Concerning regularization of pairs of vortices, D. Smets and J. Van Schaftingen [31J also 
studied the following problem 

s 2 Au=(u- q -^ lnI) p + -(g-£lnI-<, in Q, 
u — 0, on dQ, 

and obtained the exact asymptotic behavior and expansion of the least energy solution by 
similar methods for (11.81) . This method is hard to obtain multiple non- vanishing pairs of 
vortices solutions. 

In this paper, we approximate stationary vortex solutions of Euler equations (II. ip with 
multiple non-vanishing pairs of vortices solutions by stationary classical solutions. Our 
main result concerning (jl.ip is the following: 

Theorem 1.1. Suppose that Q C R 2 is a bounded simply- connected smooth domain. Let 
v n : dQ — > R be such that v n G L s (dQ) for some s > 1 satisfying j dn v n = 0. Let 
k\ > 0, k,J > 0, % = 1, • • • , m, j = 1, . . . , n. Then, for any given C 1 -stable critical point 
(xit, ■ ■ ■ , x+ x±x, ■ ■ ■ ,x~ m ) of Kirchhoff-Routh function W defined by (11.6p . there exists 
Sq > 0, such that for each e G (0,£o)> problem (II. ip has a stationary solution v e with 
outward boundary flux given by v n , such that its vorticities uf satisfying 

supp(u+) C UZi B ( x i,e, Ce) for x^ £ G Q, i = 1, • • • , m, 

supp(uj;) C U] =1 B(x~ £ , Ce) for x~ £ G Q, j = 1, • • • , n 



and as e — )■ 



REGULARIZATION OF POINT VORTICES 



5 



Remark 1.2. The simplest case, corresponding to pairs of vortices (m = n = 1) was studied 
by Smets and Van Schaftingen [31] by minimizing the corresponding energy functional. 
In their paper as e — > 0, W(x+ e ,x]~ e ) — > sup W^x^xf). Even in the case 

m = n = 1, our result extends theirs to general critical points (with additional assumption 
that the critical point is non-degenerate or stable in the sense of C 1 ). The method used in 
[31] can not be applied to deal with general critical point cases. The method used here is 
constructive and is completely different from theirs. 

Remark 1.3. In this case that m = n = 1 suppose that (x^, x^J is a strict local maxi- 
mum(or minimum) point of Kirchhoff-Routh function W(x + , x~) defined by ( Q)J) , state- 
ment of Theorem 11.11 still holds which can be proved similarly (see Remark II. 5p . Thus we 
can obtain corresponding existence result in [31] . 



Theorem 11.11 is proved via considering the following problem 

' -e^u = YZi Xn t (u-q-i In \)\ - EU Xuj fa " £ 1* J " * G 

u = 0, iG dfl, 

(1.10) 

where p > 1, q G C 2 (f2), f2 C M 2 is a bounded domain, Qf(i = 1 ■ • • , m) and f). (j = 
1 • ■ • , n) are mutually disjoint subdomains of f2 such that G and x^ G Qj . 

Theorem 1.4. Suppose q G C 2 (fi). Then for any given > 0, nj > 0, i = 1, • • • , m, j = 
1, . . . , n and for any given ^-stable critical point (x^, ■ ■ ■ , x+ # , x^, • • • , x~J of Kirchhoff- 
Routh function W defined by (11.61) . t/iere exists eo > 0, such that for eache G (0,£o); fll.lOj) 
/ias a solution u s , such that the set = {x : u £ (x) — ^ In - — g(x) > 0} CC Qf, i = 

1, • • • ,m, = {x : u £ (x) — In i — g(x) > 0} CC Qj, j = 1, • • • , n and as e — > 0, each 
shrinks to xf^ G fi. 

Remark 1.5. For the case m = n = 1, suppose that (x^, x^fj is a strict local maximum(or 
minimum) point of Kirchhoff-Routh function W(x) defined by (jl.6p . then statement of 
Theorem 11.41 still holds. This conclusion can be proved by making corresponding modifi- 
cation of the proof of Theorem II .41 in obtaining critical point of K(z) defined by (14.11) (see 
Propositions 2.3,2.5 and 2.6 in [12] for detailed arguments). 



As in [13], we prove Theorem II .41 by considering an equivalent problem of (jl.lOp instead. 
Let w = and 5 = £ (T^h)~*~, then ( Jl.lOp becomes 

-5 2 Aw = YZi Xn+ (w-k+- jtr\^ x )) P + ~ Xqj (]%t\<l( x ) ~ K J ~ W Y + > in ^ 
w = 0, on dQ 

(1.11) 
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We will use a reduction argument to prove Theorem 11.41 To this end, we need to 
construct an approximate solution for (II. lip . For the problem studied in this paper, the 
corresponding "limit" problem in IR 2 has no bounded nontrivial solution. So, we will follow 
the method in [HI EE] to construct an approximate solution. Since there are two parameters 
5, e in problem (ll.lip and two terms in nonlinearity, which causes some difficulty, we must 
take this influence into careful consideration and give delicate estimates in order to perform 
the reduction argument. For example we need to consider (sf s , ■ ■ ■ , s , s is , • • • , s~ s ) and 

( a t,si " ' a m,<5' " " " > a n,s) together in Lemm a Q 

As a final remark, we point out that problem (II. lip can be considered as a free boundary 
problem. Similar problems have been studied extensively. The reader can refer to [FT] IT3| 
[III [US EH [23] for more results on this kind of problems. 

This paper is organized as follows. In section 2, we construct the approximate solution 
for (II. lip . We will carry out a reduction argument in section 3 and the main results will be 
proved in section 4. We put some basic estimates used in sections 3 and 4 in the appendix. 

2. Approximate solutions 
In the section, we will construct approximate solutions for (II. lip . 

Let R > be a large constant, such that for any x e Q, Q CC Br(x). Consider the 
following problem: 



-5 2 Aw = (w-a) p + , mB R (0), 
w = 0, ondB R (0), 1 ' ' 

where a > is a constant. Then, (12. ip has a unique solution Ws >a , which can be written as 



^aln^/lnf, s s < \x\ < R 

where 4>{x) = <fi(\x\) is the unique solution of 

-A0 = 0f, 0>O, ^^(BjfO)) 

and ss € (0, R) satisfies 

\n(s 5 /R) 

which implies 



^| In (5|C^-i)/2 \ a J 
Moreover, by Pohozaev identity, we can get that 

/ p+1 = ^^|0'(l)| 2 and / p = 2vr|0'(l)|. 

For any z6 0, define Ws jZ ,a{ x ) = W$, a {x — z). Because Ws^ z ,a does not vanish on dVL, we 
need to make a projection. Let PWs jZ , a be the solution of 
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Then 



where g(x, z) satisfies 



-5 2 Aw = (W s , z , a - a) p + , in Q, 
w = 0, on dVl. 



PW 5)Z , a = W S , x ,a ~ t^r9(x, z), (2-3) 
In — 



—Ag = 0, in fi, 
q = In | R | , on d£l. 



It is easy to see that 



g(x, z) — In R + 2irh(x, z), 

where h(x,z) = —H(x,z). 

Let Z = ( Z+,Z ~), where Z+ = (zf,---,z+), Z+ = (zf, • • • , z~). We will construct 
solutions for (11.111) of the form 

m n 
i=l j=l 

where zf , zj G Q, > 0, ajj > for i — 1, • • • , m, j = 1, • • • , n, cjg is a perturbation 
term. To make cos as small as possible, we need to choose a^, ajj properly. 
In this paper, we always assume that zf, zj e Q satisfies 

d(z?,d£l) > g, d(zj,dtt) > g, \zf-z^\>g L , i, k = 1, ■ ■ ■ , m, i ^ k 

(2.4) 

\zj ~ h I > £ L > \ z t ~ z j I ^ Q L ' 3> l = • ■ ■ ' n ' J ^ 
where o > is a fixed small constant and L > is a fixed large constant. 

Lemma 2.1. For 5 > small, there exist (sj'^Z), ■ ■ ■ , s~g m (Z), sJ^Z), ■ ■ ■ ,sJ n (Z)) and 
( a s,i( z )i ■ • ■ j a tm( z )i a s,i( z )i ■ • • , a J,n( z )) satisfying the following system 

i»/H) ( ,+)-VHV (1) = _!__ j = i, . . . , m ( 2 .5) 

p/<P-D( s -)-^-^> (1) = — ^— , j = 1, ■ ■ • , n (2.6) 

ln{ Sj /R) 

and 

a^i 1=1 s ; 

(2.7) 
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Sj Pft k=l s k 

{21 



CL A = K- — 



where G{x, y) = In - g(x, y) for x^y. 



Since the proof is exactly the same as in Lemma 2.1 in [13] , we omit it here therefore. 

% 



To simplify our notations, for given Z = (Z+,Z n ), in this paper, we will use af v sf 
instead of af^Z), sf^Z). From now on we will always choose 



<5,1> ' " J a <5,m> a 8,V ' " J a <5,7 



and (s$ >1} ■■■ , sj m , s si , ■■■ , s s>n ) such that hold. For (aj 1? • • • , aj m , a si , ■■■ , a s J 

and (s~g V ■ ■ ■ , s^ m , sj v ■ ■ • , sj n ) chosen in such a way let us define 



P£z,i = PWs, z +,a t , P&j = PW S ^ (2.9) 
Remark 2.2. As in [13] . we have the following asymptotic expansions: 

1 1 



In 4" lu- 
st, £ 



In 


lne| 


hie 


2 



+ ° -ft— nr > i = V --.^ ( 2 - 10 ) 



dine ln^ ^ In ^ ^ In ^ V 



a^i £ 1=1 



In e\ 



Moreover, a 5j - and have similar expansions. 
To simplify notations, set 



a=l (3=1 



(2.11) 



Of^), = 0(^3), z = m,/ i = 1,2. (2.12) 



Then, we find that for x 6 -B is + (z 4 + ), where L > is any fixed constant, 
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|lne| d >*t < a ss K ' ln-f v ' * ' 1 |lne | 



(*) - ^ - j^Ktf.tf) -^({Dg(zt,zt),x- zj) + 0(\x- zt\ 2 )) 

=H W*> - K * - iff - w <zW) ' x " 



(st? 



ln-f" v 11 l ' In 4" ' W V l ln ^l 

s S,i s S,i 

and for fc ^ i and a; G B Ls + (zf), by fl272|) 



S S,k S S,k 

JtrG(4, 4) + j^b (^ + , 4)> x - 4> + o ( j ln£| 



and 



2 

W*) = j^T^, z") + ^± (DG(z+, zj), x - **) + 
So, by using (12. 7p . we obtain 

^)-^)-^ + - 27r ' /(r) 



ln-§- v 1 ' 3 ' ln-§- N -vt >~*/>~ l /'^V|lne| 



|lne| 

| In -1 ^ ' ' In- 



27r a ■ 

-- w s,z+, «+ (*) - 4». - ( D <i(4), x - 4) - ( D 9(4i 4), x - 4) 



m + n — 

k + l S tk 1=1 S S,l 

+ xeB L4 {zt). 



|lne| 



(2.13) 
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Similarly, we have 

2nq(x) 



|lne| 



=W,.- „- (x) - a 5 -. + j^L (Dq(z-),x- zj) - {Dg{zJ ,zj),x - *r) 



5,z. ,a 



n — m -| 



i^j fc=l S 4 " 



(2.14) 



We end this section by giving the following formula which can be obtained by direct 
computation and will be used in the next two sections. 



+ + (x) 



dz th 
( 1/ 



\(P+l)/2 llt \x ~ zf\, zf h ~X h ( 1 \ +N 



A|#(l)j|ln 



In -5- la; — 



±12 



|lne| 



(2.15) 



3. THE REDUCTION 



Let 



w(x) 



/0(N), 



\x\ < 1, 



1 0'(1) In |x| > I. 

Then iu G C 1 (1R 2 ). Since <f>'(l) < and In \x\ is harmonic for \x\ > 1, we see that u> satisfies 



- Aw = w p + , in M 2 . 



(3.1) 

Moreover, since w + is Lip-continuous, by the Schauder estimate, w G C 2,a for any a G 
(0,1). 

Consider the following problem: 



-Av - pwl^v = 0, v G L°°(R 2 ), 



(3.2) 
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It is easy to see that J^, i — 1,2, is a solution of (13. 2\j . Moreover, from Dancer and Yan 
|15j . we know that w is also non-degenerate, in the sense that the kernel of the operator 
Lv := -Av-pw^v, v G D 1 ' 2 ^ 2 ) is spanned by {§^ } g}. 
Let Pg~ Zi , Pfzj be the functions defined in (12. 9p . Set 



dP, 



F s . z ={u:ue L p (n), I ~-^±u = 



dP> 



s,z,j 



n 9z j h 



u = 0, 



i = !,-•• ,m, j = !,•-• ,ti, /i = 1,2 



and 



z = 1, • • • ,ra, j = 1, • • • ,n, h = 1,2}. 
For any u G L p (f2), define as follows: 

= - £ t "U (-** @0 ) - 1 1 >* (^ 

i=i /i=i \ »,/» / j=i h = i \ 

where the constants satisfy 
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dP^~ dP~ 

Since f n *f ,k Qsu = 0, J n *'f a Q$u = 0, the operator Q$ can be regarded as a projec- 

k , h IJi 

tion from L p (fl) to F$ t z- In order to show that we can solve (I3.3P and (13. 4p to obtain bf h 
and b~.^ we just need the following estimate ( by ( 12.121) and ( 12.151) ): 



ifcMl ln£ lp+l \Jl n£ lp+l 

where c > is a constant, 5 ife ^ = 1, if i = fc and h — h; otherwise, = 0. 
Similarly, 

L A (%)^ = + (n^) • (3 - 6) 

si J,h l t h 11 

where c > is a constant, = 1, if j — I and h — h; otherwise, = 0. 
Set 



L s u = -5* An - f> n + ( P + z - P s - Z -4- 

i=l ^ 11/ 



p-i 

u 



and 



^=(u£ 1 s Ki (^))U(u?=iS £ .- 

We have the following lemma. 



Z 3 



Lemma 3.1. There are constants p > and 5 > 0, such that for any 5 G (0, S ], Z 
satisfying (12 .4 p . w G -E^z Q&L$u = m \ /or some L > Zarge, i/ien 



||Q 5 L 5 u|| LP(n) > — P ° (p _ 1)2 |MU°°(n)- 
| In 5 1 p 

Proof. Set Sjy • = •. In the sequel, we will use || • || p , || • ||oo to denote || • \\u>(Q) and 
II ' IU°°(o) respectively. 
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We argue by contradiction. Suppose that there are 6^ — > 0, Z N satisfying (" 12 .4[) and 
u N e E Sn ,z n with Q Sn L 5n u n = in Q \ B s z and ||rtjv||oo = 1 such that 



1 S p 

\\Qs n L Sn u n \\ p < 



hi d) A? | p 



First, we estimate hf h A1 and b. z Ar in the following formula: 



m 2 / <9P H 

2 a S N,Z N) i 



Qs n L Sn u n = L 5n u n Kh,N ~ S N A 



i=l h=l 



dz, 



i.h 



i=i h=i V hh 



For each fixed k, multiplying (13. 7p by — ^+ Zjv,fc , noting that 



5P+ ^ fc 

k,h 



(3.7) 



we obtain 




Using ( 12TT3|) . ( 12TT4|) and Lemma [OJ we obtain 
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dP> 

, q UnLsn \ dz 



Sm ,Zfj ,k 



kji 



6 N A 



fdP+ 

I °N 



8=1 



y^jPXnf ( P 5 N ,Z N P 8 N .Z X 



2 M {x)y- 1 dpi 



Z N ,k 



\\ne N \J + dz+ h 



-V»y (p- -P+ 1 I 27Tq{x) Y~ l dP Lz N ,k 

j=l x ' k,h 



UN 



P I \W SN!ztN , at k -a SNtk 



p-i fdW t 



N> z t,N' a 5 N ,k ^ a i N ,k 



dz 



k.h 



dz 



k.h 

p-i 



'JV,< 



\\ne N \ 



P 



Y 



<9P + 

Sjy ,Zn ,k 

dz + - } 

p-1 



MAT 



dP> 



Sat .Ziv.k 



dz 



-U N 



k.h 



'N 



\\ne N \* 



Using (I3.5P and (I3.6p . we find that 



&i,iv = (41 W|) 



Similarly, 



Kh,N = 0(e 2 N \lne N \) 
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Therefore, 

m 2 / p j p+ \ n 2 / an- 

E E "t,« -«a^%^ +EE W -^a%-" j 

i=i h=i V / j=i h=i \ OZ 3\h 

-P2^2^%h,N \ W 5 N ,z+ N ,ai ■ a 8 N ,i) q+ q + 

■ ''2-2- W { W ^l N ,a^ 3 ~ a SNj ) _ 



. dz~r Oz . , 

2_i \ / 2_ 1 

m 2 p | l+ I \ / n 2 _p 1 



-o(EE f y¥ + o tE E 



,jy i 

N 



Thus, we obtain 



1 A l 1 ^ 

In on p 



For any fixed define 



Let 



L N u 



ut, N {y) = u N (s Ni y + z+ N ), u jtN (y) = u N (s Nj y + z j>N ). 



,± ^2 



2nq 



p-i 



u 



5 



p-1 



( p s N ,z N ( s %iy + z tN) - p 8 + N ,z N (4^y + 4,n) - ^ + |^|) M - 



Then 



2 5 2 ~ . 

X T^T^\\ L Nut N \\p = \\L Sn U N \\ p . 
i S N,i) 

Noting that 

5 N \ „ / 1 



4, J " Vlln^r 1 
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we find that 



Ls N u N = o 

As a result, 



°N 



^ JV I , (p-1) 2 

'lnJjvl v 



L%u± =o(l), mlf(n%), 



where = {y : s%^y + z± N eti}. 

Since ||^ jy-||oo — 1> by the regularity theory of elliptic equations, we may assume that 

ufjf^uf, in CUR 2 ). 

It is easy to see that 

E ( s tr,i) 2 ( p+ /+ , + s p - / + j_ + \ + J^!_V 1 



Similarly. 



p-1 

2 / O \ P- 1 



— >■ U>+ . 

Then, by Lemma IA.lt we find that satisfies 

—Au — pw^T u = 0. 
Now from the Proposition 3.1 in 1151 . we have 



dx\ dx 

Since 



uf = cf— +tf — . (3i 



/ A( '"£»'< )*„ = 0, 
Jn dzf_ h 



we find that 



tP-i U V „,± 



which, together with (13. 8p . gives uf ^ 0. Thus 



h 

1 i ii r i i i 
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u± N ^0, mC\B L (0)), 

for any L > 0, which implies that Un = o(l) on dB Ls ± (zf N ). 
By assumption, 

Qs N Ls N u N = 0, in Q \ B SnZn . 
On the other hand, by Lemma IA.14 for i = 1, • • • , m, j = 1, • • • , n, we have 

(^U - ^ - 4 - y^) + = o, x 6 n+ \ ^(^), 



Thus, we find that 

-Au N = 0, in O \ B SNt z N . 
However, un = on <9Q and un = o(l) on dB$ Nt z N - So we have 

^ = o(l). 

This is a contradiction. 



□ 



From Lemma [XU using Fredholm alternative, we can prove, as in [T3], the following 
result: 

Proposition 3.2. QsL$ is one to one and onto from E$.z to Fg^z- 
Now consider the equation 

QsLsu = Q s lj - Q s lJ + QsRgiu) ~ QsRj(oo), (3.9) 

where 

* = £ *t {** - - 4 - W) - E - • (»•«•) 

i=l ^ I I / + j=i 

j=l V I I / + j=1 

and 
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i=l 



P£z ~ P&.Z + U ~ K i ~~ 



2nq(x) 



In el 



p s,z P&,z K i 



2nq(x) 
I lne| 



V Ptz P&.z K i 



2nq(x] 



6,2 *6,Z | ln£ | 



P S,Z - P tz -U-Kj 



p-1 



Co' 



(3.12) 



27rg(x) 
I lne| 



27iq(x) 
I lne| 



2irq(x)\ p - 1 



Using Proposition 13.2} we can rewrite ( 13. 9 p as 



to 



G s u =: (Q s L s )- l Q s (lt - IJ + Rj(oo) - Rj(u)). 



(3.13) 



(3.14) 



The next Proposition enables us to reduce the problem of finding a solution for (11. lip 
to a finite dimensional problem. 



Proposition 3.3. There is an 5q > 0, such that for any 5 G (0, 8q\ and Z satisfying (12. 4p . 
(13. 9 P has a unique solution us G E$ z , with 



INHoo = Of 5 1 hx8\' 

Proof. It follows from Lemma IA.1I that if L is large enough, S is small then 

2irq(x) 



p+ _ p- _ K + _ 
5 > z s > z 1 | In e| 



0, x G nf\B Ls +.(zf), i = 1,... ,m 



2irq(x) 
I lne| 



0, S efi:\B is -.( Z :),i = l, 



, n. 



Let 



M = £* z n<MML< 



5|ln5|V|. 



Then M is complete under L°° norm and G$ is a map from Eg,z to -E^z- We will show 
that Gs is a contraction map from M to M. 
Step 1. Gs is a map from M to M. 

For any u G M, similar to Lemma [A. 1\ it is easy to prove that for large L > 0, 5 small 
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P, 



p- + K + 2nq ^ 
^- j ^->--->'> - |l n£ | 

P- -P+ CO K- I 2?rg(x) 



0, xenj\B Lst {z^ 



6,i 



o, ze07\s tSj _( Z 7). 



Note also that for any u G L°°(Q), 



Q s u — u in O \ 

Therefore, using Lemma fA.lt (I3.10p - fl3.13p . we find that for any u; G M, 

Q s (lt-ls) + Q s (R+(u)-Rj(u)) 
=lj-lj + Rj(uj)-Rj(uj) 
=0, in n \ B z>s . 

So, we can apply Lemma [3.11 to obtain 

\\(Q s L 5 y l (Q s (lj - lj) + Qs(Rj(co) - Rj(u))) |U 



< 



(p-1) 2 

c|ina|-?- 

2 

c) P 



IIW-i7) + Q*W(w)-i? 



10J 



Thus, for any u G M, we have 



II^MIU = \\{Q 6 L 6 )-*Q s (lt - lJ + Rj(u)-Rj(u))\\ 
< cllnSl ^ \\Qs(li - lj + Rt(u) - Rj(co) 



— 

5P 



(3.15) 



(3.16) 



It follows from fl3.3l) - fl3.6l) that the constant br*, corresponding to u G L°°(0), satisfies 



\ i,h Jn 



dP> 



6,Z,i 



<9P 



\u\ 



Since 



lj - IJ + Rt(oo) - RJ (u) =0, in Q \ B 8 , z , 
we find that the constant 6^ corresponding to it — lj + Rti 00 ) ~ Rji^) satisfies 
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<ciinirE(E/„ 

i,h \a=l JB 



B Lsf (*+) 



dP> 



6,Z,i 



dz 



i,h 



+c\\ a 6r%l£ t f 



dP> 



S,Z,j 



dz. 



j,h \/8=l JtS Ls-^ z l3> 3,h 



\lt-li + Rt(")-RJ(") 



KCe 1 -^ \nenit - lj + B+(co) - Rj(u) 



As a result, 



\\Qs(lj-lJ + Rj(u)-Rj(uj))\\ p 
<\\lt ~ IJ + Rt(u) - Rj(u)\\ p + CJ2K 



i, h 



-5 2 A 



dP+ A 



dz. 



i.h 



j,h 



9z 7,n 



<C(\\lj\\ p + \\lj\\ p + \\Rj(u)\\ p + \\Rj(u)\\ p ). 



On the other hand, from Lemma [A. II and (I2.13p . we can deduce 



8 IIP 



a ii 



Cs 



I In el v s ^' a 
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For the estimate of ||i2/(a;)||p, we have 

27!-q(x)\p 
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IKHIU = 



1=1 



\\ne\ 



+ 



P£z Ps,z 



K 



-P\Plz-Pi.z-<- 



2nq(x)y 



\\ne\ 



<C||w 



|2 

oo 



i=l 



T,x«t( p tz-Ptz-4 



2nq(x) 



p-2 



\\nS\ 



Similarly, we have 



L 5 \\p 



o 



zEIXI ' W R 6 Mil* 



2 1 V 



\\n5\ 



Thus, we obtain 



||G*M||oo< 



C\lnS 



5p 



II^IIp + II^IIp + II^ + HII p +II^HII p ) 



(p-i) 

<C|ln5| * 



<J| Inkj- 



et; 



(ln^l^+p + |ln<5| 



Thus, Gs is a map from M to M. 

Step 2. Gs is a contraction map. 

In fact, for any Ui G M, i = 1,2, we have 

GWi - GW 2 = (QsLs^QslR+M - Rj(u 2 ) - (RjM - RjM)]. 
Noting that 



R+M = RtM = 0, in n \ UZiBl s j. CO, 

0,1 



and 



2-nq(x)\p 



(3.17) 



(3.18) 



RJM = RIM = o, in n \ u] =1 b l (z~), 

0,3 J 
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we can deduce as in Step 1 that 

(p-i) 2 

HGW1-GW2II0C < C|ln<5 l " (\\RtM-Rt(u 2 )\\ P + \\RiM-RjM\\ P ) 

ip-1 / || w l||oo Halloo 



<C5| ln^l^Hwi - w 2 ||oo < r ll^i ~~ w 2||oo- 

Combining Step 1 and Step 2, we have proved that Gs is a contraction map from M to 
M. By the contraction mapping theorem, there is an unique us G M, such that a;^ = GsUJs- 
Moreover, it follows from (I3.18P that 

Halloo < 8\ ln<5| V. 

□ 



4. Proof of The main results 

In this section, we will choose Z, such that Pg~ z — Pf z +us, where ug is the map obtained 
in Proposition 13.31 is a solution of (11. lip . 
Define 




and 

K(Z) = I (P+ z - P s - Z + us ) . (4.1) 

It is well known that if Z is a critical point of K(Z), then P s + Z — P 5 ~ z + us is a solution of 
(11. lip . In the following, we will prove that K(Z) has a critical point. 

Lemma 4.1. We have 



K{Z) = I (P+ z - P, z ) + O (|^) . 
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Proof. Recall that P+ z = £™i ^ = EJ=i *W We have 

K(Z) = I{P+ Z - P 5 - z ) + 5 2 [ D{P+ Z - P 5 - z )Du s + f / |ZM 



m 1 

^ p + 1 



^ p + 1 Jf7 1 



27tq(x)^ p+1 ~ 



2iiq(x) 



p+i 



In el 



U S - Kj + 



2nq(x) 



p+i 



27rg(x) 
ThTeT 



p+i 



Using Proposition 13.31 and (I3.15p . we find 



2nq(x) 



|lne| 



P&,z Ps,z K i 



2irq(x) 
I lne| 



p+i 



5,i 



27iq(x) 



p+1 



=0 



«i) 2 IMIc 

I lnelP 



-^<5,Z P&,Z K i 



2irq(x) 



=0 



I \ne\P, 
On the other hand, 



5 2 / DP+ z Du s = / K 



<5, 2 + a+ _ a 5,i j U 5 



7/7 



i=l J B s + (z+) 

d.k 



--0 



E- 



|lne|? 
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Next, we estimate 5 2 j n \Duis\ 2 . Note that 



8=1 ^ 



US - K- 



|lne| 



+ i=i 



m 2 



+ 3=1 



n 

E K-,^ 



+EE<>. H 



i=l h=l 

Hence, by (EH-(EH, we have 



s ^Vee^-^)- 



W \Dco s \ 2 =J2 
Jn i=1 Jn+ 



|lne| 



i=i Jn 7 

m 2 



-E 

m 

+ EE"S 

i=l ft=l 
m 



P~ — P 



5,Z -US- Kj + 



v 



i./i 



lne| 

n 2 



W<5,z+a+ - a S,i 

1 0,8 



1 3 ' *,) " 



CJ5 



W<5 



i=l 



_)_ \ 'I'd- 



p- 1 / 



In el 



i=l h=l 



8z i,i J 



m 2 



llnel? 
i h=i 1 1 



71 



=0 



^E 

VIM?- 1 



p-i / s 



to 



n 2 e|6"-| 



Other terms can be estimated as above. So our assertion follows. 



Lemma 4.2. VFe have 



□ 





d 








d 


dz~r 

J,h 


3,h 



-J {Ptz- p ^z)+0{ 



| lne|P" 

e 3 

llnel^- 1 



3 = 1, " ■ ■ ,n. 



Proof. We only give the proof of the first estimate. 
First, we have 
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dK{Z) 



riPtz-pzz+us), 



dPlz dP s - z _ dug 



dz th dz th dz th 



d 



in 



dztA P& ' z Ps ' z 



dug 



£4 



P lz ~ P 8.Z 



2vrg(x) x " 
'■' " |lne| 



p s,z p s,z K k 



2iTq(x) 



X 



dPjz dP s - z 
dz th dz th 



E 



P S,Z ~ P t,Z - US - K, + 



dz th dz th 



8,Z S,Z I | ln£ j 



Since G -E^z, we have 



^4 " 



uz k,h 



dz 



jT I w <5 = 0. 

k,h 



Differentiating the above relation with respect to z^, we can deduce 



m 2 



-5 2 A- 



dR 



S,Z,t 



dug 



./» / \ 11 >>••' -".v- 



- a. 



m 2 
n 2 

/ ™ * elft+J n 2 e i 6 

\ lne\P \ lne\P 



n 2 

EE 1 ..; 



P-l f 9W 5,z+,c 



2 A dP 8,Z,/3\ dug 



-8 A 



dat n \ dug 



dz 



a,h 



dZ lh dZ lh 



P-i f dW s,z-, a ^ da Stf) \ dug 



o 



dz lh 

(— 

\\hie\P 



dz ik) 9z t 



-i 



dz li J dz th 



On the other hand, using (13.171) (for the definition of Rg(u), see (I3.12p ). we obtain 



26 



DAOMIN CAO, ZHONGYUAN LIU, AND JUNCHENG WEI 




Other teams can be estimated as above. Thus, the estimate follows. 



□ 



REGULARIZATION OF POINT VORTICES 27 



Proof of Theorem ] 1.4\ Recall that Z = Z n ). Set 

m n m 

HZ+,Z~) =^4vrV<7U + ) -^n 2 K Jq(zr) + ]T n(4) 2 g(zt, zf) 

i=l j=l i=l 

n 

+ Y n ( K j) 2 9( z 7' Z J) - Y 7lK t K tG(z+, 4) - kkJk^G(z[, zj 

3 = 1 i^k 



i=i j=i 



Zj ). 



Note that the Kirchhoff-Routh function associated to the vortex dynamics now is 

j m 1 71 

w(z+,z-)=- K t^t G ( z t,4) + 2 Y K 7 K 7 G ( z 7> z r) 

i,k=l,ij^k j,l=X,j^l 
m 1 n 

+ 2 Y^) 2 H{zi,zi) + -^( K t) 2 % )Z -) 
i=i j=i 

m n m n 

-YY K t^G{zf, Z j) + ^tM4) - Y k 7^o( 

i=l j=l i=l j=l 

Recall that h(zi, Zj) = —H{z i) Zj), it is easy to check that 

tm n 
i=l j=l 

Hence, &(Z+, %n) an d %n) possess the same critical points. 

By Lemma 14. 1| 14.21 and Proposition \A.2\ \A.3\ we have 



"7 12 



d 2 


In 


lne| 




me 


3 



+ 



d 2 


In 


lne 




lne 


3 



and 

dK{Z) 5 2 d$(Z) 
Oz± h \kis\ 2 dzf h 

Thus, the existence of a C 1 -stable critical point of Kirchhoff-Routh function W(Z) implies 
that K(Z) has a critical point. 

Thus we get a solution ws for (II. lip . Let w e = ^^-Ws, 5 = e 2 , it is not difficult 

to check that u £ has all the properties listed in Theorem 11.41 and thus the proof of Theorem 
11.41 is complete. □ 



Now we are in the position to prove Theorem 11.1 
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Proof of Theorem \l.l[ By Theorem II A\ we obtain that u £ is a solution to (jl.lOp . 
Set 

v e = (V(w e - u e = V x v £ , 

p+i 



n 1 ( *4\\YYE\\ 

n 1 / 

+ E;TT*«; 



2tt c / 2 1 

/ + 

Then (v e , P e ) forms a stationary solution for problem (11.11) . 
We now just need to verify as e — > 



j 

By direct calculations, we find that 



„ m n 

.7 = 1 .7 = 1 



| lne| 



IM P /• / _ + 27r g y ^ \\ne\v f ( 2vq 
2- 1 {2^] Qt \ WS Ki \\ne\) + (2n)Pe 2 J Q - \\\ne\ 



— K- 



LsJ JzT) 
0,1 - % ' 



J=l Ls7.(z") 



E 



f M'e 2 y a +,J J Bl(0) T 

-Ef^-E^-w 

8=1 sT. j'=l s, . 

m n 

~H K h ase^O. 
i=i i=i 

Therefore, the result follows. 
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Remark 4.3. To regularize pairs of vortices with equi-strength k, we do not need Xn+ an d 

i 

X n - , that is, we only need to consider the following problem 

-e 2 Au = (u - q - £ In A)* - (q - £ In \ - u) p + , x6fi, 
u — 0, x G d£l. 
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Appendix A. Energy expansion 



In this section we will give precise expansions of I (P$~ z — P s z ) and {Ps~z ~ P& z)i 

i,h 

which have been used in section 4. 

We always assume that zf ', z~ G Q satisfies 

d{z~l, dVL) > g, d(zj, d£l) > g, \zf — z^\>g L , i, k — 1, • • • , m, i ^ k 
\zj ~Zi\> g L , \zf - zj\ > g L } j, I = 1, ■ ■ ■ , n, j ^ I, 
where g > is a fixed small constant and L > is a fixed large constant. 
Lemma A.l. For x G i — 1, 2, • • • , m and x G Qj, j — 1, 2, • • • , m, we have 

Ps,z{ x ) ~ P sA x ) > K t + x G B 44l _ Ts}i) (z+), 

Plz (x) ~ Ptz ( x ) > «J - , x e S s - (1 _ Ta _ } (zj) , 

where T > is a large constant; while 

PtzW - Pi,zi x ) < < + t^> **nt\ PsUiHst^t), 

P 8,z( x ) - P£z( X ) < K j ~ | j^p X G \ B \ j (l+(ss j r)( Z j )' 

where a > is a small constant. 

Proof. The proof is exactly same as Lemma A.l in [13]. For reader's convenience, we give 
the proof for P^ z — Pf z here. 

Suppose that x G B + r l _ Ts + \ {zf). It follows from (I2.13P and <f>[(s) < that 
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if T > is large. On the other hand, if x G flf \ B, + )»(z^~), where a > a > is a fixed 
small constant, then 

E, , R JJ 2Trq(x) , , 



i=i 



<Ca-K+ + o(l) < 0. 
Finally, if x e B {s +^{zf) \ B 5 + ( i +T(s + y){4) for some h then 



In 



R 



"^ln ^ ~ °* + ° I hT^ 

+ Mi + r(4,f) /(»+„ 

— — + °kr-H <0 ' 

S <5,i \ s a,i 



if T > is large. Note that by the choice of a, B s + ^ 1+ ^ s + y^(z^) D 5 s +.(i + t( s +.)*)(X + ) f° r 
small 5. We therefore derive our conclusion. □ 

Proposition A. 2. We have 



^ j "ln« + 4(ln*) 2 Ll^+L^) +L |ln e ||ln5| 

£ v e ' \ i=l j = l / i= l 1 11 £ 1 

llnellln^l ^ (ln^) 2 ^ (ln^) 2 

J = l 1 11 £ 1 1 = 1 V £ / J=l V £ / 

7TP44G{4, Zj) y, n5 2 Kj^G(zr, Zj) 

tf. (ln^) 2 On-) 2 

" " 2*pKtKjG{zt,zj) / £ 2 ln|ln g | \ 

trpr ( in f) 2 v im 3 

where C is a positive constant. 
Proof. Taking advantage of ( 12. 3p . we have 
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TYL TYL 

* 2 1 nptz - p, z )\ 2 = E E J n (w s , 4 , 4k - at k y + p s % 

+ E E / (^-^, - "vY+p&j - 2 E E / - <*)+*W 

1=1 j=l JQ j=l i=l JQ 



3 

First, we estimate 



JB S + <X ) JB S + <X ) 

8,i 5,i 

In zt- Jb + fel 



. 'S s+ (Zi) 

2(p+i) p , r _ 2 P 



(4) (4/ / + 4 '"'KJ 2 / <A" 

Vs «,i' ^Bi(O) ^5,/ M(0) 

B l( o) llln^r 1 
>r(p+l)« 2 K,) 2 , 2rf 2 K,) 2 2rf 2 (a+> 2 , + +1 , „ I 



m ~r Ks s,i / J BUG) 



(ln4) 2 In 4- (In ^ ' ^ ' + U \ \ \ne^ J 



"5,, 

Next, for k ^ i, 



at 



2p 



*>« <S,fc 

~^w G ^ Zi) LJ +0 to 



>5,k' m 7f. JBl ^ 

0,1 



l n _R_||l n _^| ^ ' *' Vl ln£ l P+1 
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Moreover, we have 



(W s , 4 , 4a - <,,) p + I W St2 - a - - r^jrg{x,zj) 



z. ,a: . o,r/ 4- 

i 1 6,1 ' i 



. ,a, . i 

j 0,3 in 



i,3 



2p 



R . 

8,3 s 6,i 



\x - z- 



<■> \p-i a 8,j( s ti) 2 ^ f - +N 

R ^\ Z j 1 Z i ) 



In 



b p + 



27i5 2 at,ac ■ 

* 5 > J G(zt,zr) + 



ln-?-||ln-§- 

s S,i s S,j 



Bi(0) 

(— 



(41 

\ne\P +1 



By Lemma [O and (ESJ) , 




Other terms can be estimated as above. So, we have proved 
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7T 



4 |ln-^| 2 | In-?- 1 



ln-^| 2 



7r(p + 1) ^ 2 (%) 2 vr5 2 (a 5j ) 2 tt^(^ , z j )5\a & j f 



4 |ln-^| 2 |ln-§-| 



ln-^| 2 



+ |ln-?-||ln^-| ^ 

fe ^ 4.i s t.k l +i 



ln-?-||m-§-| 

s S,l s 6,i 



™ " 2nG(zt, zj)S*at ti ai d ^ , ^ 

2-^2-^ I It, JLII It. JLl 9 I 



m 'a+) 2 



t=l J=l 

7T5 2 



ln^||ln-§- 

S,i S.j 



Kj) 2 



i=i 



In 



H 12 



f_£L_ 



Thus, the result follows from Remark 12.21 



Proposition A. 3. We have 



□ 



d j i p-f. p_ \ 47r 2 ^/,+ ag(z+) 2tt5 2 «) 2 ^(^,4 

•| : • 1 I' <5,Z 1 8,z) - '■>■ — 



| In e\ | In— | dz. 



(ln^ 



K\2 

e ' 



i./i 



Q z + ' 2^ ' " ' 



(ln^ 



R\2 



3*7 



<5 2 In | lne| 
| In ^ | 3 



J_r( P+ P -) d J^ll i ^{KTfd 9 {zJ,zj) 

dzj- h [ 5 > z s ' z) |lne||lnf| dz~- h (lnf) 2 dzj- h 
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Proof. Direct computation yields that 
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For k i, from (12.131) . we have 




Using (I2.13P , Lemma IA.1I and Remark 12. 2\ we find that 
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since 



/ 0^ 1 (|x|)0'(|x|) r A = --|0'(l)|. 
Other terms can be estimated as above. Thus, the result follows. 

□ 
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